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Different Aspects of Relativistic Toda Chain 

S.KharchevQ, A. MironovQ, A. ZhedanovQ 



^ . ABSTRACT 
(D 

' We demonstrate that the generahzation of the relativistic Toda chain (RTC) 

is a special reduction of two-dimensional Toda Lattice hierarchy (2DTL). We 
also show that the RTC is gauge equivalent to the discrete AKNS hierarchy 
and the unitary matrix model. Relativistic Toda molecule hierarchy is also 
^ ' considered, along with the forced RTC. The simple approach to the discrete 

, RTC hierarchy based on Darboux-Backlund transformation is proposed. 

o 

(N 

1 Introduction 

Q> , Since the paper of Ruijsenaars [1], where has been proposed, the relativistic Toda chain 

(RTC) system was investigated in many papers [2]- [4]. This system can be defined by the 
equation: 

rT' 'ln= (l + eq„){l + eq„+i)-—^ 

^ . 1 + e^exp(q„+i - qn) 

n , • Ul , • ^ exp(gn-gn-i) 

' which transforms to the ordinary (non-relativistic) Toda chain (TC) in the evident limit 

^ , e — > 0. The RTC is integrable, which was discussed in different frameworks (see, for 

example, [2]- [4] and references therein). The RTC can be obtained as a limit of the general 
Ruijsenaars system [1]. 

In this paper we are going to review different Lax representations of the RTC, and to 
establish numerous realtions of it with many well-known integrable systems like AKNS, 
unitary matrix model etc. It is also shown that the RTC hierarchy can be embedded to 
the 2DTL hierarchy [5]. 

Besides, we discuss the forced RTC hierarchy and its finite analog, the relativistic Toda 
molecule. At the end of this short paper we describe the simple approach to discrete evo- 
lutions of the RTC which is based on the notion of the Darboux-Backlund transformations 
and can be considered as a natural generalization of the corresponding notion in the usual 
Toda chain theory. 

This short paper is minded as a brief review. For more details we refer the reader to 
our lengthy paper [6] . 
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2 Lax representation for RTC 

Let us describe the Lax representation for the standard RTC equation. The usual procedure 
to obtain integrable non-hnear equations consists of the two essential steps: 

i) To find appropriate spectral problem for the Baker- Akhiezer function(s). 

ii) To define the proper evolution of this function with respect to isospectral deforma- 
tions. 

Lax representation by three-term recurrent relation. In the theory of the usual 
Toda chain the first step implies the discretized version of the Schrodinger equation (see [7] , 
for example). In order to get the relativistic extension of the Toda equations, one should 
consider the following "unusual" spectral problem 

<I>„+i(2:) +a„$„(z) = z{$„(z) -h6„$„_i(z)} , n G Z (2) 

representing a particular discrete Lax operator acting on the Baker- Akhiezer function 
This is a simple three-term recurrent relation (similar to those for the Toda and 
Volterra chains) but with "unusual" spectral dependence. 

As for the second step, one should note that there exist two distinct integrable flows 
leading to the same equation (|^). As we shall see below, the spectral problem (H) can be 
naturally incorporated into the theory of two-dimensional Toda lattice (2DTL) which de- 
scribes the evolution with respect to two (infinite) sets of times {ti, t2, ■■■) , (i_i, t„2, ••■) 
(positive and negative times, in accordance with [5]). Here we describe the two particular 
flows (at the moment, we deal with them "by hands", i.e. introducing the corresponding 
Lax pairs by a guess) which lead to the RTC equations (0). The most simple evolution 
equation is that with respect to the first negative time and has the form 

^^fl..*,.-. (3) 

with some (yet unknown) i?„. 

The compatibility condition determines i?„ in terms of a„ and 6„ i?„ — ^ and leads 
to the following equations of motion: 

dttn _ K _ K+i . dbn _ ^ 1 _ ]_\ 
dt. 

In order to get (|^), we should identify 

a„ = exp(-ep„); ^« = -e^ exp(g„ - g„_i - ep„) (5) 

Note that in this parameterization the " Hamiltonian" i?„ depends only on coordinates 
q„'s: Rn = -e^exp(g„ - g„_i). 

Performing the proper rcscaling of time in (jj) we reach the RTC equation (|l]). 

As we noted already, the evolution (^), which leads to the RTC equations is not the 
unique one. The other possible choice leading to the same equations is 

^ = ~5„($„ - z$„_i) (6) 

Oil 

The compatibility condition of (||) and (^ gives the equations 

^-^ = -a„(&„+i - 6„) ^-^ = -fo,i(fe,i+i - 6,1-1 + a„_i - a„) (7) 
oti Oti 
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where 



This leads to the same RTC equation (|l|). 

2x2 matrix hax representation. The same RTC equation can be obtained from 
the matrix Lax operator depending on the spectral parameter [3] (generalizing the Lax 
operator for the TC [7]). Then the RTC arises as the compatibility condition for the 
following 2x2 matrix equations: 

L<«^,„=i^„H . ^=A„V<» (8) 

V -eexp(-g„ + ep„) J \ Vn' J 

_ f e^exp((?„ - 9„_i + ep„_i) -eC"^ exp((j„) \ ^^^^ 
V eC"^exp(-(7„_i + ep„_i) 1 - + e'^ J 
One can easily reduce these equations to the system (||) and (|^). 

To conclude this section, we remark that L-operator (^), which determines the RTC 
is not unique; moreover, it is not the simplest one. Indeed, we shall see that there exists 
the whole family of the gauge equivalent operators, which contains more "natural" ones 
and includes, in particular, the well known operator generating the AKNS hierarchy. From 
general point of view, the whole RTC hierarchy is nothing but AKNS and vice versa. 

3 RTC and unitary matrix model, AKNS, etc. 

Now we are going to describe the generalized RTC hierarchy as well as its connection 
with some other integrable systems. We start our investigation from the framework of 
orthogonal polynomials 

Unitary matrix modeL It is well-known that the partition function r„ of the 
unitary one-matrix model can be presented as a product of norms of the biorthogonal 
polynomial system [8] . Namely, let us introduce a scalar product of the form |^ 

<A,B>=<f 1^ exp| J2 (imz"" - t.mZ~n}Mz)B{-) (n) 

?ri>0 

Let us define the system of polynomials biorthogonal with respect to this scalar product 

/i„(5„fc (12) 

Then, the partition function r„ of the unitary matrix model is equal to the product of 
/i„'s: 



n-l 

Tn = 

k=0 



l[hk , To = l (13) 

A:=0 

The polynomials are normalized as follows: 

$„(z) =z" + ... + 5„_i, $*(z) = z" + ... + 5;:_i, 5_i = 5*1=1 (14) 

^The signs of positive and negative times are defined in tfiis way to get the exact correspondence with 
the times introduced in [5]. 
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These polynomials satisfy the following recurrent relations: 
$„+i(z) = z$„(z) + 5„z"$*(z-i) 



(15) 



and 



1 ~ SnS* (16) 



The above relations can be written in several equivalent forms. First, it can be presented 
in the form analogous to (||): 

<^n+i - ^<f„ = z ($„ - -^(1 ~ 5„-i5:_i)$„_i| (17) 

K+1 - -^K = z-^ [k - -1^(1 - s,,^,si_,)K-i] (18) 

From the first relation and using d) and (|), one can immediately read off 

exp(-ep„); = -e^ exp((7„ - (7„_i) (19) 



Thus, the orthogonality conditions ( p^ ) lead exactly to the spectral problem for the RTC. 
We should stress that equations (|l^) , (|l^) can be derived from the unitary matrix model. 

Using the orthogonal conditions, it is also possible to obtain the equations which de- 
scribe the time dependence of $„, $* . Differentiating (p^ with respect to times ti , t-i 
gives the evolution equations: 

(see general evolution equations with respect to higher flows below). The compatibility 
conditions give the following nonlinear evolution equations: 

dSn _ „ hn+l dSn _ „ ft-n+l , . 

Oti hn Ot-i hn 

dS* _ _c* dS* _ hn+l 

dh hn ' dt-1 ~ ""+1 hn ^^-^^ 

As a consequence, in the polynomial case, 

dhn dhn 

— ~S„Sj^_ihn', = Sn-lSj^hn (24) 

These are exactly relativistic Toda equations written in somewhat different form. Indeed, 
from (11), (|2|), (|2|) and one getsf] 

hn+l 

5 '"^ ^ " '"+0 Tife (25) 

hn 



^The same equation holds for t_i-flow. 
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On the other hand, the RTC is a particular case of the 2DTL hierarchy. Indeed, let us 
introduce the key objects in the theory of integrable systems - the r- functions, which are 
defined through the relation /i„ = T„+i/r„. Then, with the help of (p2|)-(p^, one can show 
that the functions t„ satisfy the first equation of the 2DTL: 

OtiOt_i logT„ = 2 (26) 

Therefore, it is natural to assume that the higher flows generate the whole set of non-linear 
equations of the 2DTL in spirit of [5]. This is indeed the case. 

This completes the derivation of the RTC from the unitary matrix model. 

RTC versus AKNS and "novel" hierarchies. Now let us demonstrate the corre- 
spondence between RTC and AKNS system. We have already seen that the orthogonality 
conditions naturally lead to the 2x2 formulation of the problem generated by the unitary 
matrix model: 

" \ ( Z z'^Sn \ 







A* 

^ n 




d 






\ ^ n 








$* ] 

^ n / 



ji-i 




(28) 



(29) 



(Equations (pS), ( p9| ) follow from (pO[)-(pl|) and the original spectral problem (p^). Put 
$„ = z'^/'^^'^Fn, = 2""/^+^/^^. Then the spectral problem (^ can be rewritten 
in the matrix form 

L'^^^^-^T,, = .f„+i , .F„ ^ 1^ ^ (30) 

where 

L<™=(^^^ , Q^z'" (31) 

This is the Lax operator for the discrete AKNS [9]. Obviously the evolution equations 
(^), ( p9| ) can be written in terms of F„, F* as 



dti 



A^^T^ , = ( Zr-' 1 (32) 



Note that after introducing the trivial flow 
dJ^n 




AZ'^n , AZ'=\ ^ (33) 




, 4") = I (34) 
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we get the difference non-linear Schrodinger system (DNLS) [9] (see also [7]) generated by 
the "mixed" flow 



dT-[dto 9Li ^" ^-^^ 



f + SnSj^_i (, 



C — 1 C^n 



(35) 



fndeed, from the compatibility conditions for (|31|), ( |35| ) or, equivalently, just from (22) 
(along with the trivial evolution dtoSn = '^Sn , dtoS* — ~2S* ) one gets the discrete 
version of the nonlinear Schrodinger equation: 



dS 

Qrp — ^{Sn+l — 23,1 + Sn-i) + SnS^[Sn+l + Sn~l) 



(36) 



Note also that the "novel" hierarchy of [10] is equivalent to the RTC (and, therefore, to 
the AKNS hierarchy) as well. Namely, the Lax operator in [10], i.e. 



Ln — 



Z + UnVn 



1 



Lr. 



(1) ^ 






n 






(2) 


H 




n / 







defines the recurrent relation of the form ( |17[ ) 



Un-1 



(1) 



Z (/)[ 



(37) 



(38) 



thus revealing the connection with the RTC. Comparing ( |17| ) and (38) leads to the identi- 
fication Un — Snhn , Vn — : whcrc /i„'s Satisfy ([l6|). Moreover, from ( |37| ) and 
( p!5| ) it is easy to see that 4>n^ = $„; 4>n^ — -^^2"$* — S*_i^rij and, therefore, L„ can 
be obtained from LJf^^^' by the discrete gauge transformation: 



/ 



T — TJ r(AKNS)rT-l. 



Un 



1/2-1/4 



1 






'1* 






K 


/in / 





(39) 



Evolution equations (p^)- (24) in terms of new variables u„, ti„ have the form 

dun 



dun 
'dh 
dv„ 
dti 



Unvl 



dVn 



Un-1 
1 + U„-lU„ 
Vn+1 



(40) 



dt-i 1 + UnVn+l 

and easily reproduce the usual AKNS equations in the continuum limit since 

{dto - dfi - dt_^)un = -{Un+l - 2Un + Un-l) + + U^Vn + ■ 



(41) 



{dto ~ t^tl ~ dt_-,)vn = {vn+1 - 2Vn + - {vl + vl_^_i)un + ■■■ 

We conclude with the remark that the operator LJf' in (||) is also gauge equivalent to 

£^(AKNS) (ggg [6]). 

Non-local Lax representation. There is another form of the recurrent relations 
which is non-local (i.e. contains all the functions with smaller indices) but instead expresses 
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^n{z) through themselves. This form is crucial for dealing with the RTC as a particular 
reduction of the 2DTL. Let us introduce the normalized functions 

rn{z)^<fn{z) , r*{z-')^^Ki^-') (42) 

From (|l7|), ( p^ ) one can show that in the forced and fast-decreasing cases some proper 
solutions satisfy the equations 

n ^* 

zVniz) = Vn+liz) - SnK ^ -J^Vkiz) = CnkVkiz) 

z-^V:,{z-') = ^P:+,{z-')^S: Sk^.VUz-') ^ CknPtiz-') 

iT'n , 

k— — oo 

This expression is correct for general (non-polynomial) Vn and V* provided the sums run 
over all integer k. In the polynomial case, the sums automatically run over only non- 
negative k. The last representation of the spectral problem will be useful to determine the 
general evolution of the system. Indeed, these relations manifestly describe the embedding 
of the RTC into the 2DTL [5,11], which is given essentially by two Lax operators (£ and 
£). 

RTC as reduction of 2DTL. In order to determine the whole set of the evolution 
equations, one can use different tricks. For example, one can use embedding ( ^3| ) of the 
system into the 2DTL, making use of the standard evolution of this latter [5]. Let us briefly 
describe the formalism of the 2DTL following [5] . In their framework, one introduces two 
different Baker-Akhiezer (BA) Z x Z matrices W and W. These matrices satisfy the linear 
system: 

i) the matrix version of the spectral problem: 

CW^WA , CW = WA-^ (44) 

ii) the matrix version of the evolution equations: 



-fi = (r")+w , ^ = (£™)+w 

dW — m dW — m 

''''' =(/:")_w, ^ = (/:")_w, m = i,2. 



(45) 



dt-m ' dt- 

where "L x "L matrices C and C have (by definition) the following structure: 

^ (46) 
C^Y. diag[c,(s)]A' ; c_i(s) ^ 

i>-i 

Here diag[6i(s)] denotes an infinite diagonal matrix diag(. . . 5i(— 1), 6^(0), 6i(l), . . .) ; 
A is the shift matrix with the elements A„fc = 5n.k-i and for arbitrary infinite matrix 
^ = I]i6Zdiag[ai(s)]A* we set 



{A)+ ^ diagk(s)]A' , (A)_ = ^ diag[a,(s)]A^ 

i>0 i<0 



(47) 
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i.e. (^)+ is the upper triangular part of the matrix A (including the main diagonal) while 
is strictly the lower triangular part. 

Note that (pq) can be written in components as 



Cnk = K+i,k + bk-nin)0{n - k) , C„k ^ c^iin)Sn~i,k + Ck^n{n)9{k ~ n) ; n,keZ 

(48) 

The compatibility conditions imposed on (^),(|4^) give rise to the infinite set (hierarchy) 
of nonlinear equations for the operators C, C or, equivalently, for the coefficients bmin), 
c,„(n). This is what is called 2DTL hierarchy. 

From (|43|), one gets two matrices 

Cnk=6n+l,k-'^SnS*k^^0in~k) , fc, 71 G Z (49) 

hk 

Cnk = -^Sn-i.k - Sn-iSie{k - n) , k,neZ (50) 

which have exactly the form (^). Now using the technique developed in [5], one can get 
the whole evolution of the RTC hierarchy. However, it can be also easily obtained in the 
framework of the orthogonal polynomials. 

Evolution and orthogonal polynomials. Differentiating the orthogonality condi- 
tions with respect to arbitrary times, one can obtain with the help of (^3|) the evolution of 
polynomials Vn and V* [6]: 

^ = -[{cn-uvk; ^ = [{r'^)-uvk 

[{c^hUPt; ^ = [(/:")+]fe„n (51) 



ai,„ dt 



4 Forced RTC hierarchy 

RTC-reduction of 2DTL. Let us formulate in some invariant terms what reduction of 



the 2DTL corresponds to the RTC hierarchy. Return again to the Lax representation (43) 
embedding the RTC into the 2DTL. Using (jl^), one can easily prove the following identities 

1 _ 1 ^ Y^SkSthk = K-h^^, (52) 
f-^ hk riN "Af-i , 

k—n k—n 

Because of these identities, the matrices £ and jC^ have zero modes ~ 5^-1 and Sl_i/hk 
respectively. Therefore, one could naively expect that they are not invertible and get (using 
(H)) that 

{C£)nk = Snk ? ^ " ; {CC)nk = 5nk -^-^ (53) 

" — oo hk 

Since the reduction is to be described as an invariant condition imposed on C and £, these 
formulas might serve as a starting point to describe the reduction of the 2DTL to the RTC 
hierarchy only if their r.h.s. does not depend on the dynamical variables. It seems not to 
be the case. 
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However, these formulas require some careful treatment. Indeed, the formulation of 
the 2DTL in terms of infinite-dimensional matrices, although being correct as a formal 
construction requires some accuracy if one wants to work with the genuine matrices since 
the products of the infinite-dimensional matrices should be properly defined. In fact, this 
product exists for the "band" matrices, i.e. those with only a finite number of the non-zero 
diagonals, and in some other more complicated cases (of special divergency conditions). 



One can easily see from (43) that the RTC Lax operators do not belong to this class. 
Therefore, equations (|5^) just do not make sense in this case (this is why the interpreta- 
tion of the general RTC hierarchy in invariant (say, Grassmannian) terms is a little bit 
complicated) . 

However, in the case of forced hierarchy, some of the indicated problems are removed 
since one needs to multiply only quarter-infinite matrices and, say, the product LC, always 
exists. Certainly the inverse order of the multipliers is still impossible. Therefore, only the 
first formula in ( ^3|) becomes well-defined acquiring the form 

(££)„,. = 5„fc (54) 

This formula can be already taken as a definition of the RTC-reduction of the 2DTL in 
the forced case as it does not depend on dynamical variables. Now we will show how this 
definition is reflected in different formulations of the 2DTL. 

Determinant representation. One can show (see [6]) that r- function of the forced 
hierarchy r„ = 0, n < has the following determinant representation [12] 



Tn{t) = det 



d\A-dt^y f A(z,Hexp|^(t„ 



- t_™w"™) } dzdw 



where the matrix A has the form A(z, w) — ^^^^ S{z — w ^) [6]. 

Indeed, let us rewrite the orthogonality relation ( [l^ ) in the matrix form. We define ma- 
trices D and D* with the matrix elements determined as the coefficients of the polynomials 
$„(z) and $*(z) 

M^) ^ E^^.-^"'' ^ E^^.-^"' (56) 

3 j 

Then, (|l|) looks like 

DCD*^ =H (57) 

where superscript T means transponed matrix and H denotes the diagonal matrix with 
the entries Ca = and C is the moment matrix with the matrix elements 



i J— 0,...,n— 1 

(55) 



I I ^ ^^^^^^m _ I (58) 

•^T ""^^ lm>0 J 

Let us note that D {D*'^) is the upper (lower) triangle matrix with the units on the diagonal 
(because of (|l^)). This representation is nothing but the Riemann-Hilbert problem for the 
forced hierarchy. Now taking the determinant of the both sides of (ISTI), one gets 



n-l 



det = /ifc = r„ (59) 

k=Q 

due to formula (p^. The remaining last step is to observe that 

a,^dl{~dt_,y ( ^exp|E(^™^"-i-™^"'")| (60) 
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Tn{t) = det 



n X n 



(61) 



This expression coincides with (55). One can also remark that the moment matrix Cy is 
Toeplitz matrix. This proves from the different approach that the RTC-reduction is defined 
by the Toeplitz matrices. 

5 Relativistic Toda molecule 

General properties. Let us consider further restrictions on the RTC which allows one 



to consider the both products in (53). Namely, in addition to the condition r„ 0, n < 



picking up forced hierarchy, we impose the following constraint 

r„ = 0, n> N (62) 

for some N. This system should be called iV — 1-particle relativistic Toda molecule, by anal- 
ogy with the non-relativistic case and is nothing but RTC-reduction of the two-dimensional 
Toda molecule [13,14]|. 

sl{N) Toda can be described by the kernel A{z,'w) 

N 

^(z,H = ^/W(z)gW(H (63) 

k 

where f^^\z) and g'^'^^z) are arbitrary functions. From this description, one can immedi- 
ately read off the corresponding determinant representation (^5|) . 

Indeed, equation ( p6| ) and condition ( |62| ) implies that logro and logrjv satisfy the free 
wave equation dt^dt^^ logro = dtidt_i logrA? = 0. Since the relative normalization of t„'s is 
not fixed, we are free to choose tq = 1. Then, To(i) = 1, rAr(t) = x(ii)x(^-i)i where x(ti) 
and x(t_i) are arbitrary functions. 2DTL with boundary conditions (^ was considered in 
[13]. The solution to ( p6| ) in this case is given by [14]: 

r„(t) = det dl;\-dt_,y-\i{t) (64) 

with 

N 



T^t) = Y,a^^\t)-a^''\t_,) (65) 

k=l 

where functions a^'^^t) and a*^'')(t_i) satisfy 

det9,7ia«(i) = x(t), det(-at_J'-ia«(t_i) = x{t-i) 

This result coincides with that obtained by substituting into ( p5[ ) the kernel A{z, w) of the 
form (|63|). 



Lax representation. In all our previous considerations, we dealt with infinite-dimen- 
sional matrices. Let us note that the Toda molecule can be effectively treated in terms of 
N X N matrices like the forced case could be described by the quarter-infinite matrices. 
This allows one to deal with the both identities (|5^ ) since all the products of finite matrices 
are well-defined. 



''Sometimes the Toda molecule is called non-periodic Toda [15]. It is an immediate generalization of 
the Liouville system. 
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To see this, one can just look at the recurrent relation ( p3[ ) and observe that there exists 
the finite-dimensional subsystem of {N) polynomials which is decoupled from the whole 
system. The recurrent relation for these polynomials can be considered as the finite-matrix 
Lax operator (which still does not depend on the spectral parameter, in contrast to (^). 
Indeed, from ( ^3| ) and condition (|6^), i.e. hn/hn-i — as n > N (the Toda molecule 
conditions in terms of 5- variables read as Sn — S* = 1 ior n > N — 2 or n > 0), one can 
see that 



zVNiz) = PN+iiz) -VNiz), zVn+i{z) =Vn+2(z) -Vn+i{z) etc. 



(66) 



i.e. all the polynomials Vn with n > N are trivially expressed through Vn- Therefore, 
the system can be effectively described by the dynamics of only some first polynomials 
(i.e. has really finite number of degrees of freedom). Certainly, all the same is correct 
for the star-polynomials although, in this case, it would be better to use the original 
non-singularly normalized polynomials $*. 

Now let us look at the corresponding Lax operators (|49|)-(|5C|). They are getting quite 
trivial everywhere but in the left upper corner of the size N x N [6]. Therefore, one can 
restrict himself to the system of N polynomials Vn, n — 0, 1, . . . , iV — 1 and the finite 
matrix Lax operators (of the size N x N) . 

Now one needs only to check that this finite system still has the same evolution equa- 
tions (pT|). It turns out to be the case only for the first N — 1 times. This is not so 
surprising, since, in the finite system with A'^ — 1 degrees of freedom, only first N — 1 time 
fiows are independent. Therefore, if looking at the finite matrix Lax operators, one gets 
the dependent higher fiows. On the other hand, if one embeds this finite system into the 
infinite 2DTL, one observes that the higher fiows can be no longer described inside this 
finite system. Just the finite system is often called relativistic Toda molecule. 

To simplify further notations, we introduce, instead of Sn, S*, the new variables s„ = 
(— )"+^5'„, s* = {—)'^~^^S*. Then, one can reahze that the Lax operator can be constructed 
as the product of simpler ones £ = Cn^n-i ■ ■ ■ ^i, where Ck is the unit matrix wherever 
but a 2 X 2-block: 



V 



1 



\ 



Gk = 



Sk 



Skst - 1 



(67) 



Analogously, C — Ci . . . Cn~iCn with 



V 



Gk 



1 



Gk 



Sl. 



-1 



1 - Skst Sk 



(68) 



One can trivially see that CkCk = CkCk = 1, and, therefore, one obtains CC = CC = 1 
(cf. (H). 

From these formulas, one obtains that det C = det C — 1 which reminds once more of 
the sl{N) algebra. More generally, the factorization property of the Lax operators opens 
the wide road to the group theory interpretation of the RTC molecule - see [6]. 



6 Discrete evolutions and limit to Toda chain 
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Darboux-Backlund transformations. Now we are going to discuss some discrete evo- 
lutions of the RTC given by the Darboux-Backlund transformations and their limit to the 
usual Toda chain. One can easily take the continuum limit of the formulas of this section to 
reproduce the TC as the limit of the RTC, both with the standard continuous evolutions. 

The discrete evolution equations in the RTC framework were recently introduced by 
[3] in a little bit sophisticated way. Here we outline the simple approach based on the 
notion of the Darboux-Backlund transformation (DBT). More details will be presented in 
the separate publication [16]. 

Let discrete index i denote the successive DBT's. The spectral problem now can be 
written as follows: 

$„+i(i|z) + a„(i)$„(i|z) = z|$„(i|z) -I- fe„(i)$„_i(i|z)| (69) 

Let us define the first forward DBT (treating it as a discrete analog of (||)): 

$„(z + l|z) = $„(^|z) + a,W(^)$„_l(^|z) (70) 

where an\i) are some unknown functions. One requires that $„(i-t-l) satisfies the same 
spectral problem as ( |6^ but with the shifted value of i — > i — 1. Then the compatibility 
condition gives the equations of the discrete RTC: 

an(« + 1) = an-l(«) — ; b„{l + l) ^ b,a-l[l) — TYs— - (71) 

a„_l(^) - n«) On-l(«) - "ri-l(«) 



an(i) + a,\ji(j) (72) 



where Zi are arbitrary constants. 

One can also introduce the discrete analog of (|^): 

$„(i 4- l|z) = (1 - oc^^\i))^^(i\z) + za(f)(i)$„_^(j|2) (73) 
where a„ (i) are some new unknown functions of the corresponding discrete indices. We 



refer to this evolution as to the second forward DBT. Substitution of ([73|) to (|69|) gives 
quite different system of the discrete evolution equations: 

a„(z -f 1) = a„(z) ^i^); 6„(* + 1) = fe„_l(^)^^ (74) 

\-Oi\'(i) a„-i(*) 

/X , / N 1 — cSn' [i] 1 , ^ 

an (*) + &„(*) ^y-^ = ^— - (75) 

OLV(y) 1-q;;|i(«) 

This system is the discrete counterpart of the continuum system (|^). 

Actually, in [3], four different discrete systems of the RTC equations were written. From 
our point of view, the additional evolutional systems result from the backward Darboux- 
Backlund transformations which are complimentary to those described above [6]. 

Continuum limit. Introducing some discrete shift of time A > , one can rewrite 
all the equations describing the first forward DBT as follows: 

anit + A) ^ an-i(t) — 
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6„(i + A) = &„_i(t) ^"^^^ ""'If^ (77) 

After the rescahng — *■ (7 A, one gets from (fz^ ) al^^* — gA6„/a„ and, in the continuum 
hmit, A ^ the last two equations lead directly to (Q). 

The analogous equations can be written for the second forward DBT but with Zi 
rescaled as Zi — *■ ^a- It is clear that, in the limit A ^ , one reproduces the continuum 
equations (0). 

Limit to Toda chain. Now let us make the following expansion (compare with (||)) 
a„(i) ~ 1 - epn{i) ; 6„(i) ~ -e^i?„(i) 

(78) 

z~l + eA ; Zi 1 + eXi 

Introduce also functions ^'„(i) 

$„(i) ~ e"*„(i) (79) 



It is easy to see that ([70|) leads to the forward DBT for the non-relativistic Toda chain if 
one identifies 

aW(*)~e^„(z) (80) 
Indeed, in the limit e — > 0, one gets the standard Toda spectral and evolution equations. 

There also exist some other interesting limits [6] leading to the modified discrete Toda 
equations [3]. 

Concluding remarks 

From the point of view of studying the RTC hierarchy itself, the most promising repre- 
sentation is that describing the relativistic Toda hierarchy as a particular reduction of the 
two-dimensional Toda lattice hierarchy. However, even this quite large enveloping hierar- 
chy is still insufficient. Loosly speaking, the Toda lattice is too "rigid" to reproduce both 
the continuous and discrete flows of the RTC. Therefore, one should embed the RTC into 
a more general system which admits more natural reductions. 

This is done in the forthcoming publication [16], where we show that the RTC has a 
nice interpretation if considering it as a simple reduction of the two-component KP (Toda) 
hierarchy. It turns out that, in the framework of the 2-component hierarchy, the continuous 
AKNS system, Toda chain hierarchy and the discrete AKNS (which is equivalent to the 
RTC how we proved in this paper) can be treated on equal footing. 

We acknowledge V.Fock, A.Marshakov and A.Zabrodin for the discussions. A.Z. is also 
grateful to V.Spiridonov, S.Suslov and L.Vinet for stimulating discussions. The work of 
S.K. is partially supported by grants RFFI-96-02-19085, INTAS-93-1038 and by Volkswagen 
Stiftung, that of A.M. - by grants RFFI-96-02-16347(a), INTAS-96-2058 and Volkswagen 
Stiftung. 
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